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Abstract 

The thermal free energy and pressure of dispersion interaction between two graphene sheets 
described by the Dirac model are calculated using the Lifshitz formula with reflection coefficients 
expressed via the polarization tensor. The obtained results for a pristine graphene are found to 
be in agreement with computations using Coulomb coupling between density fluctuations. For 
a graphene with nonzero mass gap parameter a qualitatively different behavior for the free en- 
ergy and pressure is obtained. The Lifshitz formula with reflection coefficients expressed via the 
polarization tensor is used as a test for different computational approaches proposed in the liter- 
ature for modeling the response function and conductivity of graphene at both zero and nonzero 
temperature. 
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I. INTRODUCTION 



The van der Waals and Casimir interactions, which are known under the generic name of 
dispersion forces,- are caused by the vacuum and thermal fluctuations of the electromagnetic 
field. At shortest separations of a few nanometers dispersion forces are usually referred 
to as the van der Waals forces. At larger separations, when the relativistic retardation 
becomes important, it is customary to speak about the Casimir forces. In the last few years 
the fluctuation induced forces attracted much attention in the literature^^"- due to their 
prospective applications in both fundamental physics and nanotechnology. Specifically, a 
lot of experiments has been performed^"- on measuring dispersion forces between metallic, 
dielectric and semiconductor surfaces spaced at separations from a few tens to a few hundreds 
nanometers. 

Recently, special attention has been directed to carbon nanostructures, such as one-atom- 
thick graphene sheets, carbon nanotubes, fuUerenes etc. which possess unique mechanical, 
electrical and optical properties.-i^ These properties appear to be particularly promising to 
provide the basis for future carbon-based nanoelectronics. Keeping in mind that elements of 
nano- and microelectromechanical devices are separated by distances of the order of tens or 
hundreds nanometers, the dispersion forces acting between them are gaining in importance. 
For this reason, a lot of papers has been devoted to calculations of the van der Waals and 
Casimir forces between two carbon nanostructures and between a carbon nanostructure and 
a regular material body.—"— Particular attention has been given also the the Casimir- Polder 
interaction of different atoms and molecules with carbon nanostructures.—"— 

The many and varied formalisms were applied to calculate the van der Waals and Casimir 
forces between two graphene sheets. Here we center our attention on the approaches con- 
sistent with the most realistic Dirac model which assumes the linear dispersion relation 
of the graphene bands at low energy^ (there is also the so-called hydrodynamic model of 
graphene^i^i^i which does not take this property into account). Specifically, we compare 
the computational results obtained using the density fluctuation approach and the ran- 
dom phase approximation,— ii^i^SiSi and using the conductivity of graphene modeled as a 
combination of Lorentz-type oscillators without account^S and with account^^ of spatial dis- 
persion. It should be emphasized that the most straightforward formalism for calculation 
of the van der Waals and Casimir forces between graphene and different substances within 
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the Dirac model is based on the Lifshitz theory and exploits the reflection coeflicients of 
the electromagnetic oscillations on graphene expressed in terms of the polarization tensor 
in (2+l)-dimensional space-time.— Using this formalism, the van der Waals and Casimir 
interactions between a graphene sheet and an ideal metal plane,— and between a graphene 
sheet and different atoms^^ and plates made of various real materials^^- were computed. The 
case of graphene-graphene interaction remained, however, unexplored within this calculation 
approach. 

The present paper is devoted to calculation of the graphene-graphene thermal van der 
Waals and Casimir interactions using the Lifshitz theory and the Dirac model for elec- 
tronic properties of graphene. The reflection coefficients of the electromagnetic oscillations 
on graphene are expressed directly through the polarization tensor without recourse to the 
concept of dielectric permittivity. All calculations are performed for both gapless pris- 
tine graphene and graphene sheets characterized by some nonzero mass gap parameter. 
Note that the Dirac-type excitations in graphene become massive under the influence of 
electron-electron interaction, substrates, defects of structure and some other effects.-i^^"— 
We calculate both the van der Waals and Casimir free energy per unit area and pressure 
as functions of separation between the graphene sheets and the van der Waals and Casimir 
pressure as a function of temperature. For a gapless graphene our results are in agreement 
with computations using Coulomb coupling between density fluctuations with subsequent 
thermal average.— The latter approach is in fact equivalents^ to the nonretarded limit of the 
Lifshitz formula with the polarization of an isolated graphene sheet described in the random 
phase approximation. It was also shown, that retardation effects are of only minor im- 
portance for graphene. For graphene sheets with nonzero mass gap parameter we obtain a 
qualitatively different behavior for the free energy and pressure, as compared to the pristine 
graphene. In this case the character of force depends on the relationship between the mass 
gap parameter and the temperature. We provide a discussion concerning the comparison 
of our results with other results obtained in the literature for graphene-graphene van der 
Waals and Casimir interactions using the Dirac model.— >^>^~— 

The paper is organized as follows. In Sec. II we begin with the Lifshitz formula containing 
the reflection coefficients derived using the Dirac model. Then we present our results for the 
free energy of graphene-graphene van der Waals and Casimir interactions per unit area and 
pressure as functions of separation and temperature. Section III contains the comparison of 
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our results with the results by others obtained at both zero and nonzero temperature. In 
Sec. IV the reader will find our conclusions and discussion. For simplicity in comparisons 
with the classical limit, we preserve the fundamental constants (the Planck constant h, the 
velocity of light c and the Boltzmann constant /cb) in all mathematical expressions. 



II. DISPERSION INTERACTION OF TWO GRAPHENE SHEETS DESCRIBED 
BY THE DIRAC MODEL 

The free energy of the van der Waals and Casimir interactions per unit area of two parallel 
graphene sheets separated by a distance a at thermal equilibrium with an environment at 
temperature T is given by the Lifshitz formulai^"- 

J^(a,T) = -^Yl / ydy {In [l-4^{^0,y)e-y]+ In [l-4^{zC,y)e-y]}. (1) 
^'^ 1=0 -^^i 

Here, C,i are the dimensionless Matsubara frequencies connected with the dimensional ones 
= 2'7ikBTl/h (where / = 0, 1, 2, . . .) by the equality (i = ^i/uc with Uc = c/(2a). The 
dimensionless variable y is connected with the magnitude of the projection of the wave vector 
on the plane of graphene, by the equality y = 2a{k\ + /c^Y^"^. The prime near the 
summation sign means that the term with / = should be taken with a factor 1/2. 

In the framework of the Dirac model at T 7^ the refiection coefficients on graphene for 
two independent polarizations of the electromagnetic field, transverse magnetic (TM) and 
transverse electric (TE), were found in Ref.— . In terms of our dimensionless variables they 
are given by2^>2^ 



rxM {iCi ,y) = -~ — — , (2) 

yUoo + 2(y2 - (f) 

iv^ - CnO^tr + 2y) - ymoo 

where the dimensionless components of the polarization tensor in (2+l)-dimensional space- 
time are expressed as IIoo,tr = Salloctr/^ through the dimensional ones and trace stands for 
the sum of spatial components and n|. 

The explicit expression for the 00-component of the polarization tensor for graphene 
with a nonzero mass gap parameter A but zero chemical potential can be written in the 
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la{y^ - Cf ) / dx 



xil — x) 



A2 + x(l-x)/(0,y) 



1/2 ^ 



+ — I dx 

F Jo 



(3) 



r 



sin(27r/x) 



X < — In h + 2cos(27r/a;)e-^'("'^"^) + g-^^^"''^"^)] - -(1 - 2x) — - , , , , , , 
27r L ^ ^ J 2^ ^cosh^(r,0,l/) + cos(27rZa;) 



A2 + C2a;(i _ 



cos{2ttIx) + e 



'9(^,(1, y) 



A2+x(l-x)/(0,y) 



cosh5f(r, C;, y) + cos(27r/x) ( 



Here, a = e^/{hc) in the fine-structure constant, A = A/{}ijiOc) is the dimensionless mass 
gap parameter, vp = vp/c ~ 1/300 is the dimensionless Fermi velocity, and r = 27iT/Tcs = 
4:nakBT/{hc), where T^s is the so-called effective temperature. The dimensionless functions 
/ anf g contained in Eq. are defined as 



g{rXi,y) = — \A' + x{l-x)f{0,y) 



(4) 



1/2 



The explicit expression for the trace of the polarization tensor is given hy—^^^^ 



Utr{tCi,y) = 8a[y^ + fiCi,y)] / dx 



X(l - X) 8" / J /r\ 

+ dx (5) 



A2 + x(l-x)/(0,z/) 



— In [1 + 2 cos(27r/x)e-5("'^''2^) + e-^3{rXi,y)^ 



X 



27T 



sin(27r/a;) 



«i^(l-2x). 

2 cosh g{T, (i, y) + cos(27r/x) 

A2 + - a;)[(l - v^fCf - v%y^] cos{2nlx) + esi^^^^'y^ 



1/2 



cosh 5f(r, Ci,y) + cos{2ttIx) 



A2 + x(l-x)/(0,y) 

The Lifshitz formula for the pressure of dispersion interaction between two parallel 
j;raphene sheets takes the form 

1=0 "^fi 

We begin with computations of the free energy of dispersion interaction between two 
^raphene sheets using Eqs. ([I])-©. In Fig. [Ha,b) we plot the computational results for the 



free energy per unit area J-" at T = 300 K normalized to the Casimir energy between two 
ideal metal planes at zero temperature, Ec{a) = —TT'^hc/{720a^), in the separation region 
(a) from 5 to 1000 nm and (b) on an enlarged scale from 5 to 100 nm. The bottom and 
top solid lines in Fig. show the ratio J^/Ec at T = 300 K for graphene sheets with the 
mass gap parameter A = 0.1 eV and OeV, respectively. In Fig. [T]^b) the quantity J^/Ec 
at T = 300 K is shown by the solid lines from bottom to top for A = 0.1 eV, 0.05 eV and 
OeV, respectively, and at T = OK by the bottom and top long-dashed lines for A = 0.1 eV 
and OeV, respectively. As can be seen in Fig. [T]^a,b), at short separations the free energy 
of graphene-graphene dispersion interaction is much smaller than the Casimir interaction 
between two ideal metal planes, but becomes relatively large with increasing separation 
distance. The nonzero gap parameter depending on its value exerts some influence on the 
free energy. From Fig. Ht^b) it is seen that already at short separations from 10 to 20 nm the 
computational results at T = 300 K differ considerably from those at T = K. This means 
that for graphene thermal effects should be taken into account not only for the Casimir 
force, but for the nonrelativistic van der Waals force as well.— 

The computations of the van der Waals and Casimir pressures between two graphene 
sheets were performed using Eqs. ([2])-([6]). In Fig. |2]^a,b) the computational results for the 
pressure P at T = 300 K normalized to the Casimir pressure between two ideal metal planes 
at zero temperature, Pc'(a) = — 7r^^c/(240a^), are plotted in the separation region (a) from 
5 to 1000 nm and (b) on an enlarged scale from 5 to 100 nm. The meaning of the solid and 
the long-dashed lines is the same as in Fig. [l]^a,b). From Fig. |2]^a,b) it can be concluded that 
the pressure of graphene-graphene dispersion interaction possesses all the same properties 
as discussed above in the case of free energy. Specifically, at T = 300 K the thermal effect 
becomes large enough at short separations of about 10-20 nm and its role quickly increases 
with increasing separation depending on the value of A. 

The top solid line in Fig. Mjo) related to the case of pristine graphene (A = 0) within 
the separation region from 5 to 100 nm is in agreement with the solid line in Fig. 2 of Ref.— 
obtained using Coulomb couplings between density fluctuations with subsequent thermal av- 
erage. As was shown, thermal effects in the pressure of dispersion interaction for graphene 
are noticeable already in the van der Waals regime, i.e., at distances of tens of nanometers 
at room temperature T = 300 K. From the comparison of the top solid and top long-dashed 
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lines in Fig. Mjo) one can calculate the relative thermal correction to the Casimir pressure 

SrPia) = (7) 

at different separations for a gapless graphene. Thus, at a = 10, 20, 50, and 100 nm one 
obtains 5tP{o) = 5.6%, 20.8%, 89.4%, and 228%, respectively [remind that for customary 
metallic and dielectric materials thermal correction becomes large only at separations of a 
few micrometers comparable with the so-called thermal length^ hc/(2kBT)]. 

Our results show that for graphene sheets with a nonzero mass gap parameter the thermal 
effect (with exception of only shortest separations) becomes even larger. Thus, from the 
comparison of the bottom solid line with the bottom long-dashed line in Fig. |2t^b) (graphene 
sheets with A = 0.1 eV) at separations a = 10, 20, 50, and 100 nm we obtain 5TP{a) = 
4.7%, 25.3%, 200%, and 893%. This makes feasible an observation of the thermal effect in 
the dispersion interaction of two graphene sheets and even an estimation of the mass gap 
parameter by the results of force measurements. 

The role of the mass gap parameter in the temperature dependence of the van der Waals 
and Casimir pressure between two graphene sheets is illustrated in Fig. [31 Here, the separa- 
tion distance is fixed at a = 30 nm and the pressure magnitudes are plotted as functions of 
temperature for A = 0.1 eV, 0.05 eV, 0.01 eV, and OeV from bottom to top, respectively. As 
can be seen in Fig. |3l for any nonzero mass gap parameter there is some temperature region 
where the pressure magnitude remains nearly constant when the temperature increases. The 
larger is the mass gap parameter, the wider is this temperature region. Thus, for A = 0.1 eV, 
0.05 eV, and 0.01 eV the pressure magnitude remains nearly constant (less than 1% increase) 
up to T = 155 K, 95 K, and 45 K, respectively. Note that the same characteristic dependence 
on A for dispersion interactions of graphene with atomic systems and material plates made 
of dielectrics and metals was found earlier.—"^ The Casimir free energy remains nearly con- 
stant if the condition ksT <C A is satisfied with a large safety margins. Under the condition 
A < ksT the thermal correction becomes relatively large. 

The results of numerical computations presented above can be supplemented by the 
asymptotic behaviors of the free energy and pressure at large and short separations (high 
and low temperatures). In the case of large separations (high temperatures) the asymp- 
totic behavior of the Casimir free energy and pressure is determined by the zero-frequency 
contribution to the Lifshitz formulas ([1]) and ([6]). Using the asymptotic expressions for the 



reflection coefficients obtained earlier,— one arrives at the following Casimir free energy 
per unit area in the large separation (high temperature) limit 



^(a,T) = - 



W(3) 



Vphc 



AaakBTln (2 cosh 



2kBT 



(8) 



where ({z) is the Riemann zeta function. Note that the second term in the square brackets 
of Eq. (IH]) is small comparing with unity due to the smallness of vp- In a similar way for the 
Casimir pressure at large separations (high temperatures) it holds 



P(a,r) 



3v%hc 



SaaksT In 1 2 cosh 



(9) 



2fcsT J J 

The asymptotic expressions ([H]) and are in a good agreement with the results of 
numerical computations shown in Figs. [D^a) and[2t^a). As an example, for A = 0, T = 300 K 
the pressure values calculated using Eq. ([9]) agree with computations in the limits of 1% at 
separations a > 370 nm. In the limits of 5% the analytic and computational results agree at 
a > 150 nm. For graphene with nonzero mass gap parameter the asymptotic expressions ([8]) 
and ([9]) become applicable starting from larger separation distances. Thus, for A = 0.1 eV, 
T = 300 K the agreement between analytic and numerical computations in the limits of 5% 
is achieved at a > 800 nm. Note that the first term on the right-hand side of Eq. (Q was 
obtained in Ref.— for a gapless graphene with A = 0. The second terms in Eqs. ([8]) and 
(|9]) provide first corrections to the previously obtained result and generalize it to the case 
of graphene with a nonzero mass gap parameter. It should be stressed that the main (first) 
terms on the right-hand side of Eqs. ([H]) and ([9]) correspond to the so-called classical limi^ 
because they do not depend on the Planck constant. 

The case of two graphene sheets interacting via the nonthermal van der Waals interaction 
is restricted to the shortest separations from 1 to 3nm. Here, using Eqs. ([1]) and ([6]) at T = 0, 
A = 0, one obtains 

E(a)= J-(a,0) = --, P(a,0) = - — , (10) 
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where the constant C is given by 

he 



C 



327r2 



oo r-y 

ydy I dC{\n 

Jo 



any 



+ ln 



22/ + aTTv/JICy), 



Numerical computations of the integrals in Eq. ( ITTj) lead to 

he 



C = 0.02101- 



2.103 X 10"^° Jm = 0.131 eV A. 



(11) 



(12) 



The result ( !T2|) is in agreement with earlier obtained^ estimation 3C ~ 0.4 eV A for the van 
der Waals pressure in Eq. ( ITOl) . In the next section the above results are compared with 
other results obtained in the literature using the Dirac model. 



III. COMPARISON OF DIFFERENT RESULTS FOR DISPERSION INTERAC- 
TION BETWEEN GRAPHENE SHEETS 

We begin with the van der Waals interaction between two graphene sheets at zero tem- 
perature. In the first paper devoted to this subjectr^ the van der Waals energy per unit 
area, as in Eq. (fTOj) . was obtained with the coefficient C equal to C = 0.288 eV A. This is 
more than twice as large as our value in Eq. ( 1121) . 

In a later paper— a smaller value for this coefficient was obtained C = 2.156 x 10^'^'^ J m = 
0.134eVA in a rather good agreement with our result f[T^ . A slightly larger value 
computed^ might be explained by slightly different value of the used Fermi velocity 
{vp = 8.73723 X lO^m/s instead of vp = lO^m/s as in our work). Although the for- 
malism used^S is nonretarded, a good agreement with our fully relativistic computations 
was achieved. This again confirms the conclusions^ that relativistic retardation does not 
play a major role for graphene. 

Now we discuss the computational results obtained at zero temperature by using the 
models for conductivity of graphene in terms of Lorentz-type oscillators.— By assuming 
that over a relatively wide range of photon frequencies up to 3 eV the graphene conductivity 
is approximately constant equal to ctq = e'^/{4h), it was found^ that the Casimir pressure is 
given by the second equality in Eq. f lTU]) . For the constant in this equality it was derived^ 
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3C = 6.88 X 10~^° Jm = 0.430 eV A leading to the following constant for the Casimir energy 
per unit area C = 2.29 x 10^^" J m = 0.143 eV A. This is a slightly larger value than was 
obtained in Eq. ( !T2|) on the basis of the Dirac model. 

The conductivity of graphene over a wider range of frequencies was also modeled-- by 
assuming that its optical properties are very similar to the in-plane optical properties of 
graphite. The latter have been mapped to a series of Lorentz oscillators with the Drude 
term over the frequency range from 0.1 to 40 eV. As a result, the Casimir pressure between 
two graphene sheets was computed using the Lifshitz formula. The computational results 
normalized to the Casimir pressure between two ideal metal planes taken from Fig. 4(b) in 
Ref.— are plotted in our Fig. HJby the gray short-dashed line. In the same figure, the gray 
solid line shows the results obtainedi^ under an assumption of constant graphene conductiv- 
ity ctq {P/Pc = 0.00529). As can be seen in Fig. HI an assumption of the frequency- dependent 
conductivity of graphene leads to significant deviations at short separation distances. For the 
sake of convenience in Fig. |l]we also present the discussed above results of Ref.— (the dotted 
line which corresponds to P/Pc = 0.00497) and our results for the gapless graphene sheets 
and for graphene with the mass gap parameter A = 0.1 eV (they are shown by the top and 
bottom long-dashed lines, respectively). For a gapless graphene it holds P/Pc = 0.00485. 

It should be emphasized that although the formalism of reflection coefficients expressed in 
terms of the polarization tensor provides a reliable test for any alternative approach, at room 
temperature the application region of the results computed at zero temperature is restricted 
to only the shortest separations below a few nanometers due to large thermal effects discussed 
in Sec. II. Because of this, below we compare our results with those computed in the literature 
taking nonzero temperature into account. 

Using the Lifshitz theory and some version of the response function of graphene in the 
random phase approximation, the Casimir free energy between two graphene sheets was 
computed^ within a wide range of separations at room temperature T = 300 K. We show 
the computational results taken from Fig. 4 of Ref.^ in our Fig. [5t^a) by the dashed line for 
the pristine graphene. For comparison purposes the solid line in Fig. |5]^a) reproduces our 
results for a pristine graphene already shown by the top solid line in Fig. [T|^a) in another 
form. As can be seen in Fig. [5]^a), the dashed line deviates significantly from the solid 
line obtained using the Dirac model and the Lifshitz theory with the refiection coefficients 
expressed via the polarization tensor (remind that the results shown by the solid line are 
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in agreement with those obtained in Ref.— ). At asymptotically large separations the free 
energy of the Casimir interaction is given by the zero-frequency term of the Lifshitz formula 
([T]) and for two pristine graphene sheets the formalism using the temperature-independent 
response function of graphene leads to J-" ~ (see Table I in Ref.^^), whereas Eq. 
demonstrates the classical limit, J-" ~ —ksT/a^, as it should be at large separations (high 
temperatures). According to Ref.— the reason for deviation from the classical limit is that 
the used formalisnt^ takes into account the direct temperature effects, as given by the finite- 
temperature Lifshitz formula, but neglects the indirect temperature effects arising from the 
temperature dependence of the dielectric response of graphene (the latter are taken into 
account by the polarization tensor at nonzero temperature). Furthermore, Ref.^^ expects 
that for graphene with high doping concentration the indirect temperature effects at room 
temperature should be negligibly small. Thus, for graphene with doping electron density 
10^^ m~^, using the formalism accounting for only the direct temperature effects it was 
found^ that J-" ~ at separations a > 1 fim in accordance to the classical limit. From this 
it was concluded that for undoped graphene sheets the formalism with neglected temperature 
dependence of the dielectric response of graphene is applicable only at vanishingly small 
separations of about a few angstroms. These suppositions of Ref.^^ concerning the origin 
of disagreement between their results and the respective results of Ref.~, coinciding with 
our results, invite further investigation. 

Now we consider computations of the thermal Casimir pressure between two graphene 
sheets with both zero and nonzero mass gap parameter.-^^ In this approach the dielectric 
properties of graphene were described via the optical conductivity calculated using the Kubo 
formalism. The computationsl results^^ for the van der Waals and Casimir pressure nor- 
malized to the Casimir pressure between two ideal metal planes are shown by the bottom 
and top dashed lines in Fig. [5]^b) for graphene sheets with A = 0.1 eV and A = OeV, re- 
spectively. These results are recalculated from Fig. 2(a) in Ref.— by using their Eq. (10). 
For comparison purposes our results for the normalized van der Waals and Casimir pressure 
are reproduced by the three solid lines from bottom to top for A = 0.1 eV, 0.05 eV, and 
OeV, respectively. As discussed in Sec. II, the top solid line is in agreement with respective 
computational results of Ref.J^ where the case of gapless graphene was also considered using 
an alternative formalism. From Fig. Mjo) it is seen that the dashed lines deviate significantly 
from the respective solid lines especially at short separations below a few tens of nanometers. 
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In the limiting case of large separation distances (high temperatures) the used formalism 
leads to one half of the result for two ideal metal planes, i.e., the obtained Casimir pres- 
sure coincides with the first term of our asymptotic expression Thus, this formalism, 
although not enough precise at moderate and short separations, satisfies the classical limit. 

IV. CONCLUSIONS AND DISCUSSION 

In the foregoing we have investigated the van der Waals and Casimir free energy and 
pressure between two graphene sheets interacting via the zero-point and thermal fluctuations 
of the electromagnetic field using the Dirac model of graphene. This was done with the help 
of the Lifshitz theory where the reflection coefficients were expressed via the components 
of the polarization tensor in (2+l)-dimensional space-time. In so doing both the pristine 
graphene and the gapped graphene were considered. 

It is common knowledge that graphene and graphene-based nanostructures are the ma- 
terials of high promise for many prospective applications in micro- and nanoelectronics and, 
more widely, in nanotechnology. Because of this, it is of high priority to have reliable theo- 
retical predictions for the van der Waals and Casimir interactions between graphene sheets 
and other carbon nanostructures spaced at separations below a micrometer. At the present 
time, suspended graphene membranes of sufficiently large area are already available.— It 
is highly probable that measurements of dispersion interaction between graphene sheets 
and other carbon-based nanostructures will be performed in the immediate future. In this 
situation the reliable and confirmed theoretical results for a simplest system, such as two 
graphene sheets, are urgently needed. 

The Lifshitz theory with reflection coefficients expressed in terms of the polarization 
tensor provides a straightforward formalism for the comparison with other approaches. We 
have calculated the free energy and pressure of dispersion interaction between two graphene 
sheets with zero mass gap parameter and arrived to the results in agreement with obtained 
earlier— using Coulomb coupling between density fluctuations. Specifically, the existence 
of large thermal effect for two graphene sheets was confirmed as well as the asymptotic 
behaviors of the free energy and pressure at short and large separations. We have also 
generalized these results to the case of graphene with a nonzero mass gap parameter. In this 
case the thermal van der Waals and Casimir interactions between two graphene sheets are 
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shown to depend on an interrelation between the temperature and the mass gap parameter. 
The results obtained were compared with some other results in the literature for graphene- 
graphene interaction at both zero and nonzero temperature. This allowed to clarify the 
regions of applicability of several approaches to the definition of response function and 
conductivity of graphene starting from the measure of agreement between these approaches 
and the Lifshitz theory with reflection coefficients found using the Dirac model. 

In the future it would be topical to investigate the van der Waals and Casimir interactions 
between graphene sheets deposited on material substrates. This subject is of much interest 
for experiments in preparation. It is of interest also to take into account deviations of the 
dispersion relation for graphene quasiparticles from linearity at high energy. 
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FIG. 1: (Color online) The normalized van der Waals and Casimir free energy of graphene-graphene 
interaction per unit area as a function of separation (a) from 5 to 1000 nm and (b) from 5 to 100 nm. 
The solid lines from bottom to top are for T = 300 K and the mass gap parameter A = 0.1 eV, 
0.05 eV, and OeV, respectively. The bottom and top long-dashed lines are plotted at T = for 
A = 0.1 eV and OeV, respectively. 
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FIG. 2: (Color online) The normalized van der Waals and Casimir pressure of graphene-graphene 
interaction as a function of separation (a) from 5 to 1000 nm and (b) from 5 to 100 nm. The solid 
lines from bottom to top are for T = 300 K and the mass gap parameter A = 0.1 eV, 0.05 eV, and 
OeV, respectively. The bottom and top long-dashed lines are plotted at T = for A = 0.1 eV and 
OeV, respectively. 
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FIG. 3: (Color online) The van der Waals and Casimir pressure for graphene-graphene interaction 
at a = 30 nm as a function of temperature. The lines from bottom to top are for the mass gap 
parameter A = 0.1 eV, 0.05 eV, 0.01 eV, and OeV, respectively. 
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FIG. 4: (Color online) The normalized van der Waals and Casimir pressure of graphene-graphene 
interaction at T = as a function of separation. The long-dashed lines show our results for A = 
0.1 eV (bottom) and OeV (top). The dotted line shows the results of Ref.— . The solid and dashed 
gray lines represent the results of Ref.^^ for a constant and frequency-dependent conductivity of 
graphene, respectively. 
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FIG. 5: (Color online) Comparison of predictions at T = 300 K for (a) the van der Waals and 
Casimir free energy per unit area where the solid and dashed lines represent the results for pristine 
graphene obtained by us and in Ref.— , respectively (see Ref.^ for possible explanations of the 
discrepancy), and (b) the normalized van der Waals and Casimir pressure where the solid lines 
represent our results for A = 0.1 eV, 0.05 eV and OeV from bottom to top, respectively, and the 
dashed lines show the results of Ref.— for A = 0.1 eV (bottom) and A = OeV (top). 
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